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^ : The recent theory of 't Hooft [ Nucl. Phys. Suppl. 68, 174 (1998)] 

Q^ I models the black hole as a system endowed with an envelope of matter 

^^ ' that obeys an equation of state in the form p ={'-/ — l)p, and acts as 

O . a source in Einstein's equations. The present paper generalizes the 't 

I [ Hooft theory so as to take into account a bulk viscosity (" in the fluid. 

brv' It is shown that even a slight positive value of (^ will suffice to yield 

complete agreement with the Hawking formula for the entropy of the 
black hole, if the value of the constant 7 takes a value that is slightly 
K> ' less than 4/3. The value 7 = 4/3 corresponds to a radiation fluid. 

1 Introduction 

The recent paper of 't Hooft [|l| is an interesting extension of the usual 
theories of Hawking radiation |0 . The 't Hooft model is static. The Hawking 
particles emitted by a black hole are treated as an envelope of matter (a static 
fluid) that obeys an equation of state, and acts as a source in Einstein's 
equations. The equation of state is conventionally written as p = (7 — l)p, 
with 7 a constant lying between 1 and 2. Here 7=1 corresponds to a 
pressure-free fluid, 7 = 4/3 corresponds to a radiation fluid, and 7 = 2 yields 
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the Zel'dovich fluid in which the velocity of sound equals the velocity of light. 
Another related work that ought to be referred to in the present context, and 
which is actually prior of that of 't Hooft, is the one of Zurek and Page 0. 

There exists evidently a very natural generalization of 't Hooft 's theory, 
namely to take into account the time dependence of the Hawking evapora- 
tion process. A black hole emitting particles necessarily has to lose mass; 
accordingly the static fluid model envisaged by 't Hooft can only be an ap- 
proximation. In the present paper we will focus attention on one particular 
aspect of the time-dependent generalization of 't Hooft's theory, namely the 
influence from a possible viscosity of the self-screening atmosphere. Although 
we do not know at present how large the viscosity of the atmospheric fluid ac- 
tually is, it seems nevertheless worthwhile to examine the physical influence 
from this factor. From a general hydrodynamic point of view it is actually 
somewhat surprising to note how little attention is usually being paid to the 
viscosity concept in relativistic or cosmological fluids; after all, we know that 
viscosity often plays an important role in fluid mechanics. In the cosmolog- 
ical context, there are a few studies of the influence from viscosity; we may 
mention the extensive treatments of Weinberg [Q, and of Gr0n 0, and there 
are a few others. 

As we will see, the incorporation of a bulk viscosity ( leads to the follow- 
ing attractive physical property of the theory: the Hawking formula for the 
entropy of the black hole, which was rather unnatural to reproduce in the 
static 't Hooft formalism, can be incorporated in the present formalism in a 
straightforward way even if the influence from viscosity is very slight. If ( is 
chosen to have some small, constant, positive value, and if the quantity /5e~" 
(see Eq.(38) below) is roughly taken to be a constant, then 7 turns out to be 
a constant, slightly less than 4/3. That is, 7 is slightly less than the value 
corresponding to a radiation fluid. 

2 Einstein's equations 

Let us thus assume the existence of a spherically symmetric atmospheric 
fluid, whose properties vary very slowly with time, in the region around the 
black hole. We recall the characteristic properties of the 't Hooft model: the 
fluid is most dense around the would-be horizon at r = 2M , but is otherwise 
present everywhere, on the outside of the horizon as well as on the inside of 



it, except at the origin where there resides a negative mass. (This negative 
mass is the price paid for the existence of the atmospheric blanket.) We write 
the hne element in the form 

ds^ = -e2-(-'*)rft2 + e^f^'-'^^dr^ + r\de^ + sin^ 9dy,^). (1) 

Einstein's equations are, with G = 1, 

Gfj,u = RfMu — T^g^iuR = SnTf^i^. (2) 

Consider first the Einstein tensor G^^, in an orthonormal basis. (Reference 
[^], for instance, gives useful expressions for this tensor.) We shall need the 
expressions for the mixed components: 

1 e-2/3 
Gl = -:^ + -^a-^P'r), (3) 

G" = -^e-("+^), (4) 

r 

r^d _ r^d _ r-fr _ (-if _ r\ /c\ 

'^Q — '^(p — '^Q — '^^ — U, \^) 

G' = -^ + ^(l + 2«'-)' (6) 

g\ = G% = -e-2"(/3 + /?2 _ ai5) + e-''\ol' + a'^ - o! 15' + ?Lz^)^ (7) 

primes and dots meaning derivatives with respect to r and t. Our basic 
assumption that the fluid, and consequently also the metric, change "slowly", 
comes into play only in Eqs. (4) and (7). The other expressions are exactly 
as in the static case. 

In accordance with our slowness assumption we can neglect the first and 
second order time derivatives of /3 in comparison to unity. Moreover, we shall 
assume that /3, when multiplied with the bulk viscosity C, stays finite. Thus, 
the time scale for our "slow" time variations is characterized by 

|/3| < 1, |/3| < 1, C/^ = finite. (8) 



We can accordingly neglect /3, /?^, and a [3 In Eq. (7). This means that 
G| = G^ reduce to the same form as for a static fluid, whereas G? = 0. 

It is of interest to relate the assumptions of Eq. (8) to the rate of energy 
dissipation in the fluid. From classical fluid dynamics it is known that the 
rate of energy dissipation per unit volume, caused by bulk viscosity, is e = 
C(V ■ u)^ (cf. Sec. 79 in |^). Thus, in our case the analogous expression per 
unit proper volume becomes (cf. Eq. (15) below): 

e = Ce^ = Cp^e-^^. (9) 

From this expression we see, however, that e is a negligible quantity in our 
approximation, because of the extra factor (3 appearing in the last term on 
the right. 

Consider next the energy- momentum tensor T^^,. As mentioned, the fluid 
acts as a source in (2). Let U^ = ([/°, W) be the four- velocity of the fluid. 
We ignore spatial derivatives, but keep time derivatives, of U^. In the static 
coordinate system we work in, we can set f/° = {—goo)'^^"^ = e~°, Uq = 
— e", [/* = 0, since the fluid is practically at rest. We introduce the scalar 
expansion 9, the projection tensor hfj_^, and the shear tensor a^^: 

^ = ^;M > V = 9,.u + U^U,, (10) 



(y^ip 



\{u,,^K + u,,^hi)-hh,, (11) 

(cf, for instance, or 0). Then, if tj is the shear viscosity and C, the bulk 
viscosity, we can write T^^ as 

T^, = pU^U, + {p-Ce)Ku-2r,a^,, (12) 

assuming constant temperature in the fluid. There are thus two viscosity 
coefficients. Usually, in cosmological applications one exploits the assumed 
isotropy of the fluid to omit the shear viscosity term. It ought to be stressed, 
however, that this point is more delicate than what is often recognized. The 
reason is that the shear viscosity under usual cosmological conditions is very 
much greater than the bulk viscosity. As an example, we mention that that in 
the universe, after termination of the plasma era at the time of recombination 
(T ~ 4000 K) one has § 

r]recomb - 6.8 X 10^ g cm'^ s"\ Crecomb - 2.6 X 10"'^ g cm'^ s"\ (13) 



The shear viscosity at this instant thus outweighs the bulk viscosity by 
about 12 orders of magnitude. Accordingly, even a very slight anisotropy in 
the fluid may easily compensate for the influence from the bulk viscosity. 

Notwithstanding this remark, we shall in the following follow cosmological 
practice and include the bulk viscosity only, thus assuming strict isotropy, in 
order to keep the formalism as simple as possible. Observing the relations 



(14) 
(15) 





r^o = «> 


rr 


= $e'^^-"\ 




U0;0 = 0, 


Ur;r 


= Pe''^-'', 


we find for the scalar expansion 








e = u[ 


+ u;^ 


= /3e-". 


We now 


define the effective pressure p: 






p = p- 


ce = p 


- CPe-^, 


whereby 


we can write the energy 


-momentum tensor as 




-^ /ij^ ^ 


pU,U, 


+ phf,u- 



(16) 



(17) 

Einstein's "energy" equation in the orthonormal frame, G^ = SnT^, leads 
to 

1 - e-2^(l - 2p'r) = Sirpr^, (18) 

which is the same equation as if the fluid were non- viscous. The analogous 
"pressure" equation, G^ = Sti'T-, leads to 

1 - e"2'^(l + 2a'r) = -8npr^, (19) 

in which the influence from viscosity is present explicitly on the right hand 
side. Finally, the equation G^ = 87rT| leads to 

e-2/3(a" + a'2 _ ^/^/ + «1^^ _ g^^^ ^20) 

when account is taken of the two first members of Eq. (8). These equations 
are the same as in the time-independent case, only with the replacement 
p —^ p. 

5 



3 Discussion on the solutions of the equations 

We insert Eq.(17) into the energy-momentum conservation equation T'^'^-^, = 
and multiply with /i^^. Since ha^U^ = we obtain 

hip,. = ~{p + p)u,,Mr (21) 

We put /i = r, and observe Ur-fi = —V^qUq = a'e° to obtain 

p=-{p + p)a'. (22) 

This is the Tolman-Oppenheimer-Volkoff equation, with the replacement 
p ^>- p. Making use of this equation, we see that Eq. (20) is a consequence of 
Eqs. (18) and (19). That is, we can henceforth consider Eqs. (18) and (19) 
to be the basic governing equations. This is the same kind of behaviour as 
in the static case. 

Let us for a moment ignore viscosity. We then see that the conventional 
state equation for cosmic fluids, 

P={l- 1)P, (23) 

when inserted into (22), yields after integration 

pg^PT _ Qonst. (24) 



In particular, if 7 = 4/3 as in the case of massless particles, we obtain 
pe"^" = const. This agrees with Eq. (2.6) in |^ (the notation is different). 

Now reinstating viscosity, we have to assume a relation between p and p 
in order to solve (22). We adopt henceforth the ansatz that most closely hes 
at hand, namely to set the deviation in pressure, C/3e~", proportional to the 
pressure itself. The proportionality constant will be called S,. Thus, we put 

P={l-i)p. (25) 

In turn, this means that we can write a natural equation of state, analogous 
to (23), also for the viscous fluid, but in terms of "tilde" variables p and 7 
instead of the usual p and 7: 

P = (7 - 1)P- (26) 



Here 

7 = 7-e(7-l)- (27) 

It is to be noted that our adoption of the hnear (barytropic) equation of 
state in the form (23) or (26) is intended to be compatible with the presence 
of the bulk viscosity. The formal similarity with the equation of state for an 
ideal (nonviscous) fluid is not to be taken as if we mean identifying the fluid 
with an isentropic fluid. As is generally known, the presence of any of the 
two possible viscosity coefficients means that there is a continous production 
of entropy in the fluid; cf. Eq. (49.6) in [Q. The state equation p = (7 — l)p 
has frequently been made use of also in earlier works on viscous cosmology; 
cf., for instance, [§,0, ^. 

Instead of (24) we now obtain, by integration of (22), 

pg^_i _ Qonst. = — . (28) 

Stt 

So far, no restriction has been made on the magnitude of ^. Let us from 
now on assume, what seems most natural, that the influence from the bulk 
viscosity is so small that we can take .^ <^ 1. It becomes then possible in 
principle to solve the governing equations (18) and (19), choosing reasonable 
values for the input parameters. One obvious possibility would be to put 
7 = 4/3, whereby 7 = 4/3 — .^/3. Let us consider this possibility first. Then 

pe^'^(l + 3e«) = ^, (29) 

87r 

to first order in S,- We can scale Einstein's equations by introducing new 
variables X and Y: 

2a 

X = — , Y = e"^. (30) 

r 

Equations (19) and (18) can now be written, to order ^, as 



X 3X2 



l-e(l + ^ln(rX)) 



rV , ^^ CY 

Y + 



l-^ein(rX) 



(31) 



(32) 



Y ' ' ' X^ 

These equations agree with Eqs. (2.9) and (2.10) in |Q in the non-viscous 
case, ^ = 0. Let the solutions in the last-mentioned case be denoted by Xq 
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and Yq- As r can now be eliminated from the equations, it is convenient to 
establish a nonlinear differential equation between Xq and Yq: 



dXn 



YM-Yo + 



CYo^ 

'xT' 



MYo-2+si: 



(33) 



This equation can be integrated in the inward direction, starting from large 
values of r for which 



.^0 - ;(i - ^: 



Fn 



'1 



2M. 



(34) 



The solution was given in [Q, and will not be further considered here. As 
mentioned, it leads to a black hole of negative mass at the origin. Let us 
instead consider the viscous case: then, the situation becomes mathemati- 
cally more complex as we cannot any longer eliminate the variable r from 
the governing equations (31) and (32). It is possible, however, to work out a 
first order perturbative expansion in ^, by writing 



X = Xo(l + eXi), Y = Yo{l + ^Y^] 



(35) 



with Xi and Yi being zeroth order quantities. From (31) and (32) we derive 
the first order "pressure" equation 



rX[ = YoYi 



CYq 
3X2 



1^1-2X1-1- -ln(rXo) 



and the first order "energy" equation 



tY; = -YoY, + 



CYq 



Fi-2Xi--ln(rXo) 



(36) 



(37) 



Once Xq and Yq are known as functions of r, we can in principle calculate Xi 
and Yi by integrating these two equations, again starting with large values 
of r and integrating in the inward direction. As initial conditions for large 
r we may set Xi = Yi = 0. If the value of the parameter ^ is known (in 
practice it has to be chosen), we can thus finally find the scaled metric from 
Eq. (35). 

However, we shall not work out the solution in this case which, as men- 
tioned, was based upon the choice 7 = 4/3, in detail. What seems to be 



of greater physical interest, is the determination of the value of ^ for which 
complete match can be obtained with Hawking's formula for the entropy of 
the black hole. In Planck units, the entropy per unit surface area is according 
to Hawking equal to 1/4. As stated by 't Hooft flj, this case requires the 
constant 7 for a non-viscous fluid to have the value 2. In other words, the 
fluid has to be of the Zel'dovich type. As mentioned in |l|], it is difficult to 
imagine ordinary matter with such a high 7 value. 

So, the central point of our paper becomes the following: for which value 
of ^, and thereby also the viscosity, do we get perfect match with the men- 
tioned Hawking formula? Actually, we do not need to make any further 
calculation to answer this question, since Eqs. (18) and (19) differ formally 
from the corresponding non- viscous equations only in the replacement p ^ p. 
Consequently, if 

P = Ip (38) 

in Eq. (19) we can carry out the same analysis as in |l| and obtain exactly 
the Hawking entropy. Since then 7 = 4/3, we obtain from Eq. (27) the 
relationship 

e = C/3e-" = ^^^^. (39) 

7-1 

Since ^ is assumed small, we can here replace a by oq. Thus, in order to 
obtain complete match with the Hawking entropy the only condition is that 
Eq. (39) has to be satisfied. The values of ^ and 7 by themselves do not 
have to be fixed. As a working hypothesis we may roughly assume that the 
quantity /3e~"° is constant; then the bulk viscosity (, as well as the parameter 
7, turn out to be constants. The constancy of 7 should be expected, in view 
of the equation of state for the fiuid, adopted as it is in the form (23). 

One may ask: is the expression (39) physically reasonable? The answer 
turns out to be in the affirmative, due to the following reason: in a realistic 
fiuid, we expext that the value of 7 in the state equation p = (7 — l)p lies 
somewhere between 1 (pressure-less fiuid) and 4/3 (radiation fiuid). This 
means that the right hand side of (39) is negative. Moreover, due to the 
emission of matter from the black hole we must have M < 0, and this 
corresponds to /? < for a given value of r. [ At large values of r, this 
follows from the last member of Eq. (34). Also in the main region of the 
atmospheric blanket, centered around r = 2M, we have /? < 0; cf. Fig. 2 in 
[0]. The only exception may be in the vicinity of the origin.] Consequently, 
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Eq. (39) yields ( > 0. The bulk viscosity turns out to be positive, as it 
should be according to ordinary thermodynamics |^ . The positiveness of the 
viscosity coefficients comes from the general thermodynamic property that 
the entropy change for an irreversible process in a closed physical system is 
always positive. 

Consider finally the question how to determine the location of the horizon 
of the time-dependent black hole. We assume that the static (i. e., the non- 
viscous) problem has been solved, so that Xq and Yq are known functions of r. 
Also, the constant C, appearing in Eq. (28), is then known. The "horizon" in 
the presence of an envelope of matter is not a singular boundary on which the 
metric diverges, but is naturally defined as the surface where the function e"^^ 
has a maximum (cf. the metric in Eq. (1)). That is, the horizon corresponds 
in the static case to the equation Yq = 0. It is given explicitly as a dashed 
line in Fig. 2 in [|I| . Mathematically, we can write the condition for the static 
horizon as 

l-Yo + ^ = (40) 



X, 



(cf. Eq. (2.10) in ll, or Eq. (32) above with ^ = 0). Moving on to the 
time-dependent, viscous case, the horizon is analogously determined by the 
equation F' = 0, or 

Y^ + aYoY^y = 0; (41) 

cf. Eq. (35). The position of the horizon is seen to be slightly displaced 
relative to the static case because of the constant small factor ^ = C/3e~°°. 
For practical purposes it is more convenient to express the condition (41) 
in such a form that the functions Xq, Yq, Xi, Yi occur, but not their 
derivatives: taking the logarithmic derivative of Eq. (35), d\nY/{dlnr) ~ 
dhiYo/{d\nr) + C^dYi/{d\nr) we get, using Eq. (37), 



CY 







Yo + ^-^Yo 



XI 



U-^-'-^-\ 



Fi-^(ri-2Xi--ln(rXo) 



0. (42) 



If the set of equations (36) and (37) has been solved, so that Xi and Yi 
are known functions of r, Eq. (42) thus fixes the position of the horizon. In 
practice, numerical work becomes necessary. 
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4 Concluding remarks 

Of course, our theory presented above is somewhat speculative. Already the 
static 't Hooft model ll| is to some extent speculative. The 't Hooft model 
is however an interesting alternative to the standard picture of a black hole, 
since it explicitly takes into account the influence upon the metric from the 
atmospheric blanket consisting of the emitted Hawking particles. Although it 
is not evident in advance that the introduction of viscosity coefficients in the 
blanket is needed in the model, this idea does not seem unreasonable to us 
in view of the general importance of viscosities in ordinary hydrodynamics. 

One important observation in our paper is that Einstein's "pressure" 
equation (19) takes the same form as the corresponding equation in the non- 
viscous case, only with the substitution p ^ p, where p is defined by Eq.(16). 
This makes it possible to use 't Hooft 's results and write down the formula 
(39), giving agreement with Hawking's entropy, directly. Thus only a slight 
bulk viscosity is sufficient to give the Hawking entropy, if the value of 7 is 
adjusted accordingly, according to Eq. (39). A small viscosity corresponds 
to a value of 7 being slightly less than 4/3, meaning that the velocity of 
sound becomes slightly less than c/v^. This seems physically reasonable. 
Moreover, we obtain a positive bulk viscosity, which is in agreement with 
ordinary thermodynamics. 

Of course, one may wonder why we find it so desirable to maintain the 
connection with Hawking's formula for the entropy of a black hole. We do 
not enter into a detailed discussion of this point, but restrict ourselves to 
mentioning that Hawking's formula seems to be the outcome of investiga- 
tions along different routes, and therefore ought to be regarded with some 
confidence in the difficult field of gravitational thermodynamics. 

Finally, it is to be noted that quantum corrections to the thermodynamics 
of the 't Hooft black hole model have recently been calculated by Nojiri and 
Odintsov [|l^. The inclusion of quantum corrections does not change the 
qualitative properties of this model. For example, the area law is found to 
be the same as without quantum corrections. 

Acknowledgement. 1 wish to thank Professor Sergei Odintsov for valu- 
able information about this problem. 
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